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The simple hydrodynamic theory of penetration by shaped charge jets is extended to include nonuniform 
jet velocity distributions and stand-off distance. The penetration of a jet with an initial linear velocity 
distribution is compared to an ideal upper limit and is found to be close to the maximum obtainable in 
practice. The conclusions concern idealized jets, without regard to the performance limitations which 


break-up and waver impose on real jets. 


I, INTRODUCTION 


HE development of a high velocity jet from a 

conical metal shell which is projected toward the 
axis by an explosive charge is illustrated by Fig. 1. 
Such jets, commonly call shaped charge jets, are easily 
able to penetrate several cone diameters into metal 
or rock targets and hence have many important 
applications. 

Figure 2 shows a jet emerging from a shaped charge. 
The decrease in the diameter of the jet as it travels is 
due to stretching caused by the velocity gradient along 
the jet. As indicated elsewhere,!? formation of the jet 
and subsequent penetration of a target are essentially 
independent processes. Thus, they may be considered 
separately ; this greatly simplifies the analysis. 

The basic theory of penetration by shaped charge 
jets was given over a decade ago.'? Using simple con- 
cepts of ideal incompressible hydrodynamics, it was 
shown that the penetration of a continuous jet of 
uniform velocity should be proportional to the length of 
the jet and the square root of the ratio of jet density to 
target density. Eichelberger,? Singh,* and Cook® showed 
that the simple hydrodynamic theory is valid for jet 
velocities greater than 0.4 cm/ysec. At lower velocities 
the strength of the jet and target become important 
and actual penetration is less than predicted.54§ 
Modifications of the simple theory have been advanced 
to account for jet and target strength at low jet 
velocities.*-§ 

With typical shaped charges, jet velocity decreases 
monotonically from front to rear. Thus the jet stretches 
as it travels toward the target and arrives with increased 
length, thereby giving greater penetration than a jet of 
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uniform velocity. As the distance from charge to target 
(stand-off distance) is increased, stretching due to non- 
uniform jet velocity is enhanced. However, deviations 
from perfect axial alignment (waver) and break-up’ of 
the jet from excessive stretching (see Figs. 3 and 4) 
cause penetration to decrease at stand-off distances 
greater than a few cone diameters. The original investi- 
gators'? recognized the effects of nonuniform jet 
velocity and stand-off on penetration and treated them 
qualitatively. In the present paper the simple hydro- 
dynamic theory is extended to include these effects. 
Since waver and break-up are not included, it is an 
investigation of what could be attained, or approached, 
if these effects could be eliminated. As in the original 
theory, the basic idealization is that the jet and target 
behave as ideal incompressible fluids at the extremely 
high pressures generated by impact of the jet. 


Il. THEORY 


Figure 5 shows a jet approaching a target. The jet at 
this instant, =0, is of length / and has a velocity 
distribution »(x). Taking £ as the initial « coordinate of 
the element which will arrive at the target (bottom of 
the hole) at time #, and assuming that the velocity »(£) 
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7 These effects are discussed elsewhere; see references 1-5. 


Fic. 1. Development 
of a shaped charge jet 
from a collapsing conical 
shell. 
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Fic. 2. Double flash x-ray photograph of a jet emerging from 
a shaped charge. The interval between the flashes was about 
MSEC. 


of the element is unchanged until impact, the displace- 
ment of the element at time / is 


w(e)=EO+P), (1) 


where P(t) is penetration. If the penetration rate is 
U (2), the velocity of an element of the jet with respect 
to the plane of impact (bottom of the hole) is 
v(£)— U (2), as indicated in Fig. 5(b). The pressure at 
the stagnation point in Fig. 5(b) is the same for the jet 
and target. Hence, assuming steady flow conditions, 
Bernoulli’s law gives 


pv—UP=pil”, (2) 


where p is the density of the jet and py is the density of 


Frc. 3. Double flash x-ray 
photograph showing jet waver. 
The interval between flashes 
was about 15 usec. 
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the target. Since v— U and U are positive (2) gives 


v= U[1+ (0/1)*]. (3) 

Differentiating (1) with respect to time, with 

dP/dt=U, and combining the resulting equation with 
(3) to eliminate U, we obtain 

di/dé+a(v'/v)t=a/», (4) 


where v’=dv/dé and a=1+ (p/p1)?. Regarding » as a 
known function of € and using v® as an integrating 
factor, (4) gives 

S41 


[o"’]}esu—[0]es=a i: ve(edz. (5) 


Ss 


Putting [vjes=v0, Lvlesp:=01, and noting that 


Fic. 4. Flash x-ray photograph 
illustrating jet break-up. 


ee 
[t]e-s=S/v0, (5) may be combined with (1) to give 
p=[ tian { dea bestov/nye—11, © 
s 


where P is total penetration and B=a—1= (p/p1)'. To 
show the effect of jet velocity distribution it is con- 
venient to substitute »= (—S)// in (6), with O<n<1 
corresponding to S<&<S-+-/. Thus we obtain 


Pail (+1/00) | oahdn—1 |+SCCon/o)*—1) (7) 


For a jet of uniform velocity v(n)=v=2, and (7) 


reduces to 
P=I(p/p1)', (8) 


which is the result obtained earlier.1:? 


PENETRATION BY SHAPE CHARGED JETS 


Ill. DISCUSSION 


The effect of jet velocity distribution on penetration 
is readily seen from (7). For fixed v1, %, and S, penetra- 
tion is maximized by giving »(n) the greatest possible 
values for 0<4<1. With the restriction dv/dn<0, this 
is accomplished by giving the entire jet, except the 
element at the rear, the velocity vo. Then (7) yields 


P=1L(B+1) (v0/1)8—1J+-SL(v0/01)?—1], (9) 


where P; denotes “ideal” penetration. The penetration 
given by (9) can never be realized in practice because 
it is for a jet with the velocity gradient concentrated at 
the rear; thus the rear of the jet is of infinitesimal 
thickness. Nevertheless, the ideal penetration is useful 
as an upper bound. 

The effect of stand-off on penetration is evident in 
(7). For a jet of uniform velocity penetration is inde- 
pendent of stand-off, and (7) (with v=21) conforms to 
this. For >, penetration increases with stand-off at 
a rate depending upon 8. As mentioned above, the 
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Fic. 5. Notation: (a) configuration at /=0, target stationary; 
(b) configuration at time ¢, bottom of hole stationary. 


increase of actual penetration with stand-off is limited 
by jet break-up and waver. 


IV. PENETRATION CALCULATION FOR 
A TYPICAL SHAPED CHARGE 


As an illustrative example the foregoing theory may 
be used to calculate the penetration of the charge of 
Fig. 6, the velocity distribution for which is given by 
Eichelberger.* Figure 7 shows the velocity distribution 
as a function of axial distance from the top of the cone. 
The upper ends of the curves correspond to the front 
of the jet and the lower ends to the rear. The curve for 
t=0 is given by Eichelberger® as the basic description 
of the jet velocity distribution; from it the velocity 
distribution at any later time can be found. The curves 
for =0 and ¢=10.2 are artificial since the point corre- 
sponding to the front of the jet is plotted to the left of 
2=0 (the top of the cone). They represent extrapola- 
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I'ic. 6. Carnegie Institute of Technology (CIT) 
standard shaped charge; see reference 3. 


tions of the experimental velocity distributions to times 
prior to the beginning of jet formation. This procedure 
is valid as long as the jet is allowed to form before it 
arrives at the target. 

The jet starts to form between z=10.2 wsec and 
20.4 usec, since the point representing the front of the 
jet in Fig. 7 passes the top of the cone in this interval. 
The straight line divides the jet approximately into 
formed and unformed parts It may be imagined that as 
a section of the velocity curve passes the line, the corre- 
sponding element of the jet forms. The entire jet is 
formed by /=40.7 usec, at which time the front of the 
jet is 20 cm from the original location of the top of the 
cone. 

The charge of Fig. 6 performs well at a stand-off 
distance of 10 cm. At this stand-off the velocity curve 
at the instant of impact with the target would lie be- 
tween the curves for /=30.5 and t=40.7 usec. Thus 
part of the jet is unformed when the front starts 
penetrating the target. 

As indicated in Fig. 7, » varies from zero at the front 
of the jet to unity at the rear. When the velocity distri- 
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Fic. 7. Velocity distribution for the CIT standard charge. The 
dashed curves represent jet elements before formation, and the 
solid curves represent jet elements after formation. 
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t= 30.5 usec 


t= 40.7 psec 


JET VELOCITY —— cm/psec 


Fic. 8. Velocity distribution for the CIT standard charge 
as a function of 7. 


bution is plotted as a function of 7, the curves of Fig. 8 
are obtained. For ‘=O and #=10.2 psec, v(m) is similar 
in shape to v(z). However, for t= 20.4, v(m) has a sharp 
break due to the vertical tangent in v(z). The same form 
is found for :+30.5 usec. Subsequently, the velocity 
distribution becomes smooth again as shown by the 
curve for t= 40.7 usec. 

The velocity curve for t=20.4 usec in Fig. 7 shows 
the front of the jet at the base of the cone. It is con- 
venient to use this particular velocity distribution in 
calculating penetration since the conventional interpre- 
tation of stand-off distance can then be retained. From 
(7), with B=1, 1=0.19 cm/sec, v9=0.75 cm/ysec, 
1=2.8 cm,® and S o'v(q)dn=0.44, we find 


P=10.24+38S. (10) 


For S=10 cm, (10) gives P=40.2 cm. This is sub- 
stantially greater than the penetration of 12 cm 
reported? for the CIT charge. The discrepancy is 
attributed to jet and target strength effects, jet break- 
up, and waver. If only the part of the jet with velocity 
greater than 0.4 cm/usec is used in the calculation, the 
penetration from (7) is 


P=3.0+0.88S. (11) 


For S=10 cm, (11) yields P=12 cm, the same as the 
observed value. But since the other part of the jet, with 
velocity less than 0.4 cm/usec, produces some penetra- 
tion, the 12-cm penetration calculated for the forward 
part of the jet must be regarded as somewhat high. The 
apparent lower penetration obtained experimentally for 


8The velocity curve shows two velocities for a given axial 
coordinate. The extrapolated fictitious jet is folded over and is 
like a string in the shape of a flattened C, which is in process of 
straightening out. The length / is the complete length, given by 
the sum of the horizontal spans of the upper and lower parts of 
the curve in Fig. 7. 


AND J. N. GOODIER 
the forward part of the jet is attributed® to jet break-up 
and waver. 


V. VARIATION IN JET VELOCITY DISTRIBUTION 
WITH TIME 


For jets of nonuniform velocity the velocity distribu- 
tion varies as the jet travels toward the target. It is 
readily shown that the variations always tend to 
produce a jet with a uniform velocity gradient. Let x 
denote the distance from the front of the jet to a typical 
element and let / be the length of the jet. At a time ¢ 
later the distance from the front to the element origin- 
ally at x is 


a =a+[oo—v(x) ]t (12) 
and the new length of the jet is 
V=Il+ (vo— me. (13) 


Dividing (12) by (13), putting n=2/l, and 7/=2'/1, 
yields 
nt[vo—2(q) }t/1 
| ar meM arn (14) 
i+ (v—»,)i/1 


The velocity gradient may be written 


dv dv dy 
—=— —. (15) 
dn’ dn dr! 
With dn/dy’ from (14), (15) becomes 
dv doy 1+(m—2)t/1 
= | i (16) 
dn’ dnL1— (dv/dn) (t/1) 
As t increases indefinitely, (16) gives 
dv 
i — (v9— 1). (17) 


Thus the velocity distribution becomes linear at 
sufficiently large ¢. 

The change in velocity gradient for smaller values of 
{may be seen by comparing the coefficient of // in the 
numerator and denominator of the bracketed factor 
in (16). The coefficient of ¢/1 in the numerator is the 
average velocity gradient. That in the denominator is 
the negative of the velocity gradient at 7 and is always 
positive if we take du/dy<0. Hence the bracketed factor 
is greater than unity for |dv/dy|<(vo—v;) and is less 
than unity for [dv/dy|>(vo—0). Thus the velocity 
distribution is always changing towards the linear form. 
This is consistent with the observation! that shaped 
charge jets generally have nearly linear velocity 
distributions. The velocity curves of Figs. 7 and 8 
illustrate this effect. (Only the portion of the jet which 
has formed should be considered.) 


PENETRATION BY SHAPE CHARGED JETS 


VI. PENETRATION OF A JET WITH A LINEAR 
VELOCITY DISTRIBUTION 


We now compare the ideal penetration of (9) with 
the penetration for a jet with an initial linear velocity 
distribution. With the linear velocity distribution 


0(m) = vo{ 1a (01/00) — 1} (18) 
and = 1, (7) gives 
P=1(v9/21)+S[ (0/01) — 1). (19) 
The ideal penetration from (9), with 6=1, is 
P.=Il (200/01) —1]+-S[ (0/01) — 1]. (20) 


For v/%,=3, a reasonable choice, the ratio of the 
penetrations from (19) and (20) is 


P/Pi= (31+2S)/ (5I+2S). (21) 


For the charge used in the example above the stand-off 
distance (10 cm) is about four times the length of the 
jet (2.8 cm) at t= 20.4 usec (see Fig. 7). If S=4l, (18) 
gives P/P;=0.85. Thus a jet with a linear initial 
velocity distribution gives about 85% of the ideal 
penetration. This indicates that there is little to be 
gained by modifying jet velocity distributions beyond 
the linear form. A more profitable approach would be 
to increase v9/?;. Since the lower limit on 0; is fixed by 
the strength of the jet and target, this must be ac- 
complished by increasing 7. 
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Vil. EFFECT OF PROJECTILE VELOCITY 


The effect of increasing the average velocity of a jet, 
perhaps by launching the charge with a rocket, is 
readily seen from (7). Adding a uniform velocity to 
v(y) increases the value of the integral in (7). However, 
the denominator 2,6 shows a greater proportional in- 
crease, since 2; is the least value of v(n). Hence the 
factor multiplying / actually decreases. The same is true 
for the factor multiplying S. Thus, for a jet and target 
which satisfy the assumptions of the theory, the addi- 
tion of a uniform velocity would cause a decrease in 
penetration. 

If, however, a portion of the rear of the jet has a 
velocity which is less than that required to overcome 
the strength of the jet and target, v, in (7) would remain 
constant with addition of a uniform velocity (until all 
of the jet is traveling at a velocity sufficient to overcome 
jet and target strength) and v» would increase. Thus, 
up to a limit, penetration would increase with addition 
of a uniform velocity. 

All these conclusions of course concern an idealized 
jet free from break-up and waver, effects which are 
important in determining the performance of real jets. 
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